function get_pertConsistency(ExpGroup)
adata_dir = '\\argon.fmi.ch\keller.g\AData\';
fb_arr=zeros(13,length(ExpGroup));
pb_arr=zeros(13,length(ExpGroup));
ids=[-6:1:7];
ids(ids==0)=[];
for fnd=1:length(ExpGroup)
    display(['Exp: ' num2str(ExpGroup(fnd)) ' *****************************']);
    [adata_file,mouse_id,userID]=find_adata_file(ExpGroup(fnd),adata_dir);
    orig=load([adata_dir userID '\' mouse_id '\' adata_file],'nbr_frames','aux_files');
    
    if length(orig.nbr_frames)>5
        order=[];
        for xnd=1:length(orig.nbr_frames)
            if orig.nbr_frames(xnd)==20000
                if xnd==length(orig.nbr_frames)
                    order(xnd)=3;
                else
                    order(xnd)=1;
                end
            else
                order(xnd)=2;
            end
        end
    else
        order=[1 2 1 2 3];
    end
    a=1:length(order);
    order(max(a(order==2))+1:end)=3;
    
    cnt=1;
    cnt2=1;
    fb=[];
    pb=[];
    for vnd=1:sum(order<3)
        if order(vnd)==1
            fb(cnt)=(str2num(orig.aux_files{vnd}(findstr(orig.aux_files{vnd},'T')+1:findstr(orig.aux_files{vnd},'.')-1)));
            cnt=cnt+1;
        else
            pb(cnt2)=(str2num(orig.aux_files{vnd}(findstr(orig.aux_files{vnd},'T')+1:findstr(orig.aux_files{vnd},'.')-1)));
            cnt2=cnt2+1;
        end
    end
    
    fb = fb_check(fb);
    pb = pb_check(pb,userID,mouse_id);
    
    for lnd=1:length(ids)
        if sum(ismember(fb(:,1),ids(lnd)))
            fb_arr(lnd,fnd)=fb(ismember(fb(:,1),ids(lnd)),2);
        end
        if sum(ismember(fb(:,1),ids(lnd)))
            pb_arr(lnd,fnd)=pb(ismember(fb(:,1),ids(lnd)),2);
        end
    end
end

curYmax=max([max(fb_arr) max(pb_arr)]);

sequ=[1 2 3 8 9 10 4 5 6 11 12 13 7];
figure;
for ind=1:length(ids)
    subplot(2,7,sequ(ind))
    hold on
    plot(fb_arr(ind,:),'r','LineWidth',2)
    hold on
    plot(pb_arr(ind,:),'black','LineWidth',2)
    axis tight
    ylim([0.8 curYmax])
    xlim([0 length(ExpGroup)])
    set(gca,'XTick',[0:1:length(ExpGroup)])
    set(gca,'XTickLabel',[0:1:length(ExpGroup)])
    drawnow
end

figure;
plot(mean(fb_arr),'r','LineWidth',2)
hold on
plot(mean(pb_arr),'black','LineWidth',2)
legend('feedback','playback')
ylim([0.8 curYmax])
xlim([0 length(ExpGroup)])
set(gca,'XTick',[0:1:length(ExpGroup)])
set(gca,'XTickLabel',[0:1:length(ExpGroup)])
drawnow